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We find that the definition of the heavy-quark fragmentation function given by Jaffe and Randall 
differs by a factor of the longitudinal-momentum fraction z from the standard Collins-Soper defini- 
tion. Once this factor is taken into account, the explicit calculation of Braaten et al. is found to be 
in agreement with the general analysis of Jaffe and Randall. We also examine the model of Peterson 
et al. for heavy-quark fragmentation and find that the quoted values of the width and of the value 
of z at the maximum are in error. The corrected values are in agreement with the analysis of Jaffe 
and Randall. 



PACS numbers: 12.39.Hg, 13.87.Fh, 14.40. Lb, 14.40. Nd 



I. INTRODUCTION 

There does not yet exist a complete first-principles cal- 
culation of the nonperturbative transition of a parton 
(emerging from a high energy scattering process) to a 
hadron. Nevertheless, a rigorous formal description of 
fragmentation functions in terms of operator matrix ele- 
ments has been available for some time Q . Additionally, 
fragmentation functions are an integral part of present 
day phenomenology. 

When the parton involved in the fragmentation is a 
relatively heavy charm or beauty quark, there are sim- 
plifications in the treatment of fragmentation functions 
that allow one to make additional theoretical progress. It 
is the case of heavy-quark fragmentation that we consider 
here. 

A theoretical description of fragmentation that is based 
on the operator-matrix-element definition of the fragmen- 
tation function in heavy-quark effective theory (HQET) 
Q has been presented by Jaffe and Randall j|. They 
show that the heavy-quark expansion of the fragmen- 
tation function must take a specific form that depends 
only on certain combinations of variables. In this pa- 
per, we point out that the Jaffe-Randall definition of 
the fragmentation function differs from the standard 
Collins-Soper definition by a factor of the longitudinal- 
momentum fraction z and that this difference is crucial 
to the proper interpretation of the Randall- Jaffe form. 

The work of Braaten et al. M models the heavy-quark 
fragmentation function with a fixed-order perturbative 
calculation in the context of HQET. The results of this 
calculation appear, at first sight, to be incompatible with 
the Jaffe-Randall form. However, after accounting for 
the difference between the Jaffe-Randall definition of the 
fragmentation function and the Collins-Soper definition, 
which is used by Braaten et al., we find agreement be- 



tween the explicit calculation and the Jaffe-Randall form. 

Perhaps the most widely-used parameterization for 
fragmentation functions is that of the model of Peterson 
et al. ||, which is based on old-fashioned perturbation 
theory. The Peterson parameterization has been fit to 
charm and beauty production data Q and subsequently 
used as input in making predictions for hadron-, photo-, 
and electro-production of charm and beauty 0. 

It has been suggested that the Peterson et al. param- 
eterization is incompatible the Jaffe-Randall form. In 
particular, it has been reported that the width of the 
fragmentation function is different in the two approaches 
H . We find that the width and value of z corresponding 
to the maximum are reported incorrectly in the original 
Peterson et al. paper |^|. The corrected values are com- 
patible with the Jaffe-Randall form. 

In Sec. II, we give an overview of the general HQET 
analysis of Jaffe and Randall. We discuss the compati- 
bility of the explicit calculation of Braaten et al. with the 
Jaffe-Randall form in Sec. III. Sec. IV contains a descrip- 
tion of the model of Peterson et al. and a discussion of 
its relation to the Randall- Jaffe form. Our findings are 
summarized in Sec. V. 



II. THE HQET ANALYSIS OF JAFFE AND 
RANDALL 

We now discuss the work of Jaffe and Randall [pj, 
which provides a QCD-based interpretation of heavy- 
quark fragmentation in terms of the heavy-quark mass 
expansion. We begin with the standard Collins-Soper 
Jy definition of the fragmentation function for a heavy 
quark into a heavy hadron: 
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where the trace is over color and Dirac indices, iV c is 
the number of colors, h{x) is the heavy-quark field at 
space-time position x, P is the four-momentum of the 
heavy hadron, n is defined by n 2 = and n ■ P — 1, and 
z = n ■ P/n ■ k, where k is the heavy-quark momentum. 
The state \H'(P)) consists of the heavy hadron plus any 
number of additional hadrons. The matrix element is 
understood to be evaluated in the light-cone gauge n ■ 
A = 0. Eqn. (|l]) is valid in d dimensions. However, in 
discussing the Jaffe- Randall analysis, we specialize to the 
case of four dimensions. 

The definition (|l|) is normalized such that f(z, /i 2 ) 
times the parton-level cross section reproduces the full 
QCD cross section, differential in z, in the collinear limit. 
By "collinear limit" we mean the limit in which Vfc 2 
and the transverse component (relative to P) of the 3- 
momentum of k are neglected in comparison to the lon- 
gitudinal component of the 3-momentum of fc.jj] This 
normalization ensures that /(z,/x 2 ) is a probability dis- 
tribution in z (Ref. jD). 



* The normalization on the right side of Eq. (|l|) can be 
understood as follows. For simplicity, we consider a frame 
in which P — (P + , P 2 /P + , Ox) has no transverse compo- 
nent, k — (k + ,k~ ,k±), with P + = zk + . In this frame, 
n = (0, l/P + ,0x) and A = P + x~, where x is the space- 
time separation of the quark and antiquark fields. In the 
collinear limit, k~ and k±_ are neglected in comparison with 
k + . In that limit, the full QCD cross section is reproduced 
by the parton-level cross section, times the Fourier trans- 
form in brackets in Eq. (|l]), times a normalization factor. 
This normalization factor includes 1/(2P + ) for the normal- 
ization of the hadron state, 2k + for the normalization of 
the quark state in the parton-level cross section, 1/4 from 
the Fierz re- arrangement I aa 'Ig/3> — * (1/4)7^3 7~ //3 / that fac- 
tors the Dirac indices, 1/N C from the Fierz re-arrangement 
Iaa'IbV — * (1 / 'N c )l a bl a ' v that factors the color indices, and 
l/k + for the quark sum-over-spinors projector (7 • k w 7~fc + ) 
in the parton-level cross section. The full QCD cross section 
is differential in the longitudinal momentum of the final-state 
hadron. Since dP + = ~k + dz, we obtain an additional fac- 
tor of k+ on conversion to a cross section that is differential 
in z. In the full QCD cross section, there is also an inte- 
gral over the transverse momentum of the hadron. This in- 
tegration corresponds, in the frame in which the hadron has 
zero transverse momentum, to an integration with respect to 
d d ~ 2 (zkj_) (Ref. |I|]). The conversion to an integration with 
respect to d d ~ 2 k± yields a factor z d ~ 2 . (The integration over 
all k± corresponds to taking the quark and antiquark fields 
at transverse separation x± — 0.) Altogether, we obtain a 
factor jz d_3 /(4iV c ), which is just the coefficient of the Fourier 



transform in Eq. (N 



We note that the definition ([!]) contains a factor z (in 
four dimensions) relative to the definition of the frag- 
mentation function used by Jaffe and Randall This 
factor z will be important in comparing with the work of 
Braaten et al. below. 

The analysis of Eq. (|l|) in HQET proceeds as follows. 
First, following the standard method for obtaining the 
heavy-quark-mass expansion, one decomposes the field 
h{x) into the sum of a large component h v (x) and a small 
component h v (x) : 



h v (x) 



— imnv-x 



P + h(x) 



h v (x) = e-" n Q v - x P-h{x), 



(2) 
(3) 



with rriQ the heavy-quark mass, P± = (1 ± ^)/2, and 
v the hadron's four-velocity. The leading term in the 
mass expansion of f(x, /x 2 ) is contained in the large-large 
combination of fields: 



-irriQ An-v 



x (0| P + h{Xn) \H\P)) (H'(P)\ P + h(0) |0) 
+ .... (4) 

Here, and throughout this paper, we use the ellipsis to 
denote terms of higher order in the hadron-mass expan- 
sion. 

In Ref. ||, it is argued that the matrix element in 
Eq. ([I|) is a dimensionless function ^(Xd). This function 
may be written in terms of its Fourier transform: 
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and ma is the hadron mass. Inserting Eq. (|J) into Eq. 
(0) , one can evaluate the integral, with the result 



f(z^) = -ai 



(7) 



A more complete analysis in Ref. || also yields the next- 
to-leading term in the hadron-mass expansion: 



/(z,^i 2 ) = z 
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My) + b(y) 



(8) 



where y = [1/z — ttlq /uih) / 5- The analysis in Ref. ||] 
does not yield a precise prediction for the functional form 
of a and b, but some general properties may be deduced. 
The function a describes, in the limit of infinite heavy- 
quark mass, the effects of binding in the heavy hadron 
on the heavy-quark momentum distribution. For a free 
heavy quark, a(y) would be a (5-function at y = 1. In a 
heavy hadron, the binding smears the heavy-quark mo- 
mentum distribution. It can be shown || that the dis- 
tribution has a width 
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IV. PETERSON FRAGMENTATION 



and a maximum at z = z max , where 

1 - z max ~ 6. (9b) 

III. A PERTURB ATIVE MODEL 

Braaten et al. Q present a QCD-inspired model for the 
fragmentation of a heavy quark into an S— wave light- 
heavy meson. In this model, the fragmentation function 
is computed in perturbative QCD (at Born level) in an 
expansion in inverse powers of the heavy-quark mass. 

Braaten et al. define the fragmentation function as the 
collinear limit of the ratio of the cross section for produc- 
ing a hadron to the cross section for producing a quark. 
As we explained in Sec. II, this definition is equivalent 
to the Collins-Soper definition (|l]) and leads to a frag- 
mentation function that is a probability distribution in 
z. 

For the projection of the Qq state onto the meson, 
Braaten et al. take the standard nonrelativistic-bound- 
state expression. For example, in the case of a 1 Sq meson, 
they assume the Feynman rule for the QqH vertex to be 
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(10) 



where -R(O) is the radial wave function at the origin. 
From the terms of leading order in the heavy-quark mass 
expansion, they obtain, in the case of a So meson, 
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(11) 



where N = 2a 2 s \R(0)\ 2 /{81nm 3 q ). 

At first glance, this result may seem to contradict the 
Jaffe- Randall analysis, which shows that the terms of 
leading order in the heavy-quark mass expansion give 
a contribution that is contained entirely in the function 
a(y) in Eq. (^J). However, the factor z in the definition of 
the fragmentation function ([!]) is crucial here. From the 
definitions of y and 8, we have z = 1/[1 — 6(1 — y)], and, 
so, we can re- write Eq. (O) as 



t . _ N (l- y) 2 fo 2 
o y 

which is of the form of a(y) in Eq. (0). 



(12) 



Finally, we examine the model of Peterson et al. 0] 
for the fragmentation of a fast-moving heavy quark Q 
with mass tuq into a heavy hadron H (consisting of Qq) 
with mass mu and a light quark q with mass m q . The 
basic assumption in this model is that the amplitude 
for the fragmentation is proportional to 1/(AE), where 
AE = Eh + E q — Eq is the energy denominator for the 
process in old-fashioned perturbation theory. It follows 
that the probability for the transition Q — * H + q is pro- 
portional to \/(AE) 2 . Taking the heavy quark's momen- 
tum to define the longitudinal axis, one can express AE 
in terms of the magnitude of the heavy quark's momen- 
tum Pq, the fraction z of the heavy quark's momentum 
that is carried by the heavy hadron, and the transverse 
momentum p± of the heavy hadron or the light quark: 



AE 



H+P 2 ± + z 2 Ph + Jm 2 q +p 2 ± + (l- z) 2 P 2 q 
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In the last line, we have set m# « wlq and neglected p\ 
relative to trig and used the definition 



e = {ml +p\)/ml 



(14) 



Multiplying 1/(AE) 2 by a factor 1/z for the longitudinal 
phase space, one arrives at the following ansatz for the 
fragmentation function || 



N 



z[l-l/z-e/(l-z)] 2 ' 
where the normalization N is fixed by the condition 
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and the sum extends over all hadrons that contain Q. 

Contrary to the claims in Ref. ||, we find that Dq(z) 
has a maximum at z = z max , with z max w 1 - ^/e, 
and a width of order y/e. Specifically we find, that 
the distance between the inflection points of Dq(z) is 
[(8 - 2V3)/3] 1/2 Ve, and the full width at half maximum 
is 2y/e. In these values of 2 m ax and the width, we have 
neglected terms of higher order in e. 

We can compare the Peterson et al. form with the 
Randall- Jaffe results by making use of Eqs. (^]) and Q), 
to obtain 



(17) 
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It has been believed that the shape of the Peterson et al. 
fragmentation function is incompatible with results ob- 
tained from heavy quark effective theory (HQET). How- 
ever, it follows immediately from Eq. (|17| ) that our results 
for the width and z max of the Peterson form are compat- 
ible with the Jaffe- Randall constraints 
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V. SUMMARY 



We have examined the Randall- Jaffe analysis of the 
heavy-quark fragmentation function in HQET and found 
that the Randall- Jaffe definition of the heavy-quark frag- 
mentation function differs from the standard Collins- 
Soper definition (Q) by a factor z. This factor is crucial 
to the interpretation of the fragmentation function as a 
probability distribution in z. We have found that, once 
this factor has been taken into account, the explicit, per- 
turbative model calculation of Braaten et al. is in agree- 
ment with the Randall- Jaffe analysis. 

We have also examined the model of Peterson et al. for 
the heavy-quark fragmentation function. Our results for 
the width and the value of z corresponding to the maxi- 
mum of the Peterson fragmentation function differ from 
the values stated in the Peterson et al. paper §. Our 
values are consistent with constraints from the general 
analysis of Randall and Jaffe. 



ACKNOWLEDGMENTS 



We thank R. Jaffe for confirming the relationship of the 
Randall- Jaffe definition of the fragmentation function to 
the Collins-Soper definition. This work was supported by 
the United States Department of Energy, High Energy 
Physics Division, under contract W-31-109-Eng-38. 



[1] J. C. Collins and D. E. Soper, Nucl. Phys. B194, 445 
(1982). 

[2] For reviews see B. Grinstein, Ann. Rev. Nucl. Part. Sci. 

42, 101 (1992); M. Neubert, Phys. Rept. 245, 259 (1994); 

A. F. Falk in Proceedings of the XXIVth SLAC Summer 

Institute on Particle Physics, edited by J. Chan, L. De- 

Porcel and L. Dixon ( NTIS, Springfield, VA, 1997), p. 43. 
[3] R. L. Jaffe and L. Randall, Nucl. Phys. B412, 79 (1994). 
[4] E. Braaten, K. Cheung, S. Fleming, and T. C. Yuan, Phys. 

Rev. D 51, 4819 (1995). 
[5] C. Peterson, D. Schlatter, I. Schmitt, and P. M. Zerwas, 

Phys. Rev. D 27, 105 (1983). 
[6] P. Nason and C. Oleari, Nucl. Phys. B565, 245 (2000); 

M. Cacciari and M. Greco, Phys. Rev. D 55, 7134 (1997). 
[7] S. Frixione, M. L. Mangano, P. Nason, and G. Ridolfi, 

Nucl. Phys. B431, 453 (1994); J. Binnewies, B. A. Kniehl, 



4 



